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Abstract 

Heterochromatin preferentially assembles at repetitive DNA elements, playing roles in transcriptional silencing, 
recombination suppression, and chromosome segregation. The RNAi machinery is required for heterochromatin assembly 
in a diverse range of organisms. In fission yeast, RNA splicing factors are also required for pericentric heterochromatin 
assembly, and a prevailing model is that splicing factors provide a platform for siRNA generation independently of their 
splicing activity. Here, by screening the fission yeast deletion library, we discovered four novel splicing factors that are 
required for pericentric heterochromatin assembly. Sequencing total cellular RNAs from the strongest of these mutants, 
cwfl4A, showed intron retention in mRNAs of several RNAi factors. Moreover, introducing cDNA versions of RNAi factors 
significantly restored pericentric heterochromatin in splicing mutants. We also found that mutations of splicing factors 
resulted in defective telomeric heterochromatin assembly and mis-splicing the mRNA of shelterin component Tpzl, and 
that replacement of fpz7 + with its cDNA partially rescued heterochromatin defects at telomeres in splicing mutants. Thus, 
proper splicing of RNAi and shelterin factors contributes to heterochromatin assembly at pericentric regions and telomeres. 
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Introduction 

Eukaryotic genomic DNA associates with histone and non- 
histone proteins to form chromatin, which is necessary for the 
spatial and temporal organization of chromosomes. A distinction is 
commonly drawn between two types of chromatin: euchromatin 
and heterochromatin. Euchromatin is less condensed and often 
associated with genes that are actively transcribed. Heterochroma- 
tin is highly condensed and often forms over repetitive DNA 
elements such as transposons. The formation of heterochromatin 
prevents expression of transposons, improper recombination of 
repetitive genomic loci, and missegregation of chromosomes during 
mitosis and meiosis, thus maintaining genome stability [1,2]. 

Histones within heterochromatin regions are usually hypo- 
acetylated and methylated at histone H3 lysine 9, which serves as a 
binding site for heterochromatin protein 1 (HP1) [3-5]. HP1 
subsequently recruits diverse proteins to regulate cellular processes 
such as transcriptional silencing, recombination suppression, and 
chromosome segregation [1]. The mechanism that attracts histone 
methyltransferases and deacetylases to repetitive DNA elements is 
under intensive study. In certain cases, sequence-specific DNA 
binding proteins are directly involved in recruitment of these 
enzymes [6-9] . Alternatively, the repetitive nature may itself be 
sufficient to trigger heterochromatin assembly [10,11]. Recent 



work has shown that DNA repeats are transiently transcribed, and 
the transcripts are processed by the RNA interference (RNAi) 
machinery into small interfering RNAs (siRNAs), which help 
target histone-modifying enzymes to repeat regions. However, the 
mechanistic details of RNAi-mediated heterochromatin assembly 
is not yet well understood [12-14]. 

The mechanisms by which heterochromatin is assembled and 
regulated have been extensively studied in the fission yeast S. pombe 
because it shares basic pathways of heterochromatin assembly with 
mammals, yet has the key advantages of facile genetics and single 
representative genes for most key families of mammalian 
chromatin-modifying factors. In this organism, heterochromatin 
is present mainly at pericentric regions, subtelomeres, and the 
silent mating-type region, all of which contain similar repeat 
sequences composed of dg and dh repeats [1]. These repeats are 
transcribed during S-phase of the cell cycle by RNA polymerase II. 
The transcripts are sliced by Ago 1 and then reverse-transcribed by 
the RNA-directed RNA polymerase complex (RDRC: Rdpl, 
Cidl2, and Hrrl) into double stranded RNAs, which are 
processed by Deri into small interfering RNAs (siRNAs). These 
siRNAs are loaded onto Argonaute siRNA chaperone complex 
(ARC: Arbl, Arb2, and Agol) and then onto the RNA-induced 
transcriptional silencing complex (RITS: Agol, Chpl, and Tas3). 
RITS is targeted to the repeat regions through base pairing 
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Author Summary 

Heterochromatin formation at specific genomic regions is 
critical for processes as diverse as gene expression and 
chromosome segregation. The formation of silent hetero- 
chromatin at repetitive DNA elements requires processing 
of transcripts by the RNA interference machinery. Curious- 
ly, factors involved in proper RNA splicing are required for 
heterochromatin assembly, and it was proposed that 
splicing factors provide a platform for the recruitment of 
RNAi complexes independently of their role in regulating 
splicing. In this study, we found several novel splicing 
factors involved in heterochromatin assembly. Our analysis 
of genome-wide splicing patterns by RNA sequencing 
showed that the mRNAs of RNAi factors are very sensitive 
to perturbations of RNA splicing machinery. Moreover, we 
found that splicing factors are critical for the production of 
a telomere shelterin component and proper telomeric 
heterochromatin assembly. Most importantly, we showed 
that introducing the cDNA versions of RNAi and shelterin 
components alleviates heterochromatin defects associated 
with splicing factor mutations. Thus splicing factors are 
required for heterochromatin assembly mainly by regulat- 
ing the proper splicing of heterochromatin assembly 
factors. 



between siRNAs and the nascent transcripts and recruits the 
H3K9 methyltransferase complex CLRC, which contains SET 
domain protein Clr4 as its catalytic subunit. H3K9 methylation 
recruits HP1 family proteins Swi6 and Chp2, which in turn recruit 
histone deacetylases (HDACs) such as SHREC to further compact 
chromatin (see review [12-14]). 

Surprisingly, in addition to these complexes, splicing factors are 
required for heterochromatin assembly at pericentric regions [15— 
17]. In fission yeast, 43% of genes contain introns [18], indicating 
the prevalence of splicing in this organism. The spliceosome and 
the splicing reactions of fission yeast are also highly conserved with 
those of higher eukaryotes, which utilize snRNAs Ul, U2, U4, U5, 
and U6, as well as over one hundred protein components (see 
review [19]). The process of splicing starts when Ul and U2 
snRNPs bind to the 5' splice site and branch point on a pre- 
mRNA, respectively. The U4/U6.U5 snRNP and the Prpl9 
complex (also known as NineTeen Complex, or NTC) are then 
recruited to form the precatalytic spliceosome. After the release of 
Ul and U4, the spliceosome is activated and the 5' splice site is 
cleaved and fused to the branch point to form a lariat structure. 
Next, 3' intron cleavage is coupled to exon ligation in a post- 
spliceosomal complex of U2, U5, and U6. The mature mRNA is 
then released and the snRNPs recycled. Most notably, tempera- 
ture-sensitive mutants prplO-1 (component of U2) and cuflO-1 
(component of U5) exhibit silencing defects at pericentric regions 
[15]. Like dcrlA, these mutants lose most siRNAs derived from 
pericentric repeats [15]. Additionally, the spliceosome associates 
with RDRC component Cidl2 [15,20], indicating a possible direct 
role of splicing factors in connecting the nascent transcripts and 
the RNAi machinery during heterochromatin formation. It was 
hypothesized that splicing factors act in the RNAi pathway 
independendy of RNA splicing because silencing defects are 
obvious when the splicing of a control tbpl mRNA is intact and 
introducing cDNAs of agol + or hrrl + was unable to rescue 
pericentric heterochromatin silencing defects of prplO-1 cells 
[15,17]. However, the possibility that splicing factors regulate 
the proper processing of RNAi factors has not been rigorously 
tested. Notably, a number of recently identified factors involved in 



RNAi (arbl + , arb2* , ersl + , and dshl + ) contain introns and might 
require splicing factors for their proper expression [21-24]. 

In this study, we performed a screen of the S. pombe nonessential 
gene deletion strain library and discovered four new putative 
splicing factors involved in pericentric heterochromatin assembly. 
We demonstrated that the phenotype of the strongest of these, 
cwfl4A, is similar to those of RNAi mutants in regulating 
pericentric heterochromatin assembly. RNA-seq analyses further 
found that cw/14A resulted in mis-splicing of a subgroup of genes, 
including a number of RNAi factors. Moreover, we showed that 
introducing the cDNAs of three RNAi factors, agol + , arb2 + , and 
ml , significantly alleviated silencing defects associated with 
cwfl4A. Furthermore, we found that the mRNA of telomere 
shelterin protein Tpzl, which is involved in telomeric silencing, 
was also mis-spliced in splicing mutants and that introducing tpzl + 
cDNA partially rescues telomeric silencing defects of splicing 
factor mutants. Thus splicing factors are involved in heterochro- 
matin assembly mainly through regulating the proper splicing of 
heterochromatin assembly factors. 

Results 

A high-throughput screen for mutants affecting 
pericentric heterochromatin 

To comprehensively identify factors required for pericentric 
heterochromatin assembly, we performed a screen of the fission 
yeast haploid deletion library for mutants that affect silencing of a 
reporter inserted into pericentric heterochromatin, otr::ade6 + 
(Figure 1A-C) [25]. In wild- type cells, otr::ade6* is silenced, causing 
red colony color on low adenine medium. However, when 
heterochromatic silencing is lost at the pericentric region, otr::ade6* 
is expressed, and colonies are white. Strains with intermediate 
silencing defects show variable degrees of pink/ red color, allowing 
rough phenotypic quantification (Figure 1C). In order to eliminate 
strains that have inherent metabolic defects causing lighter-than- 
red color, we performed a control screen with an ade6-M210 query 
strain lacking the otr::ade6 + reporter (Figure IB), which allowed us 
to filter false positives out of the screen. Our finalized list of hits is 
shown in Figure 1 C . Each hit colony was assigned a score between 
1 and 4, with 4 indicating the strongest silencing defects. 

Among the mutants identified, 25 were previously known to be 
required for pericentric heterochromatin assembly, validating the 
effectiveness of our screen. These were mutants in the complexes 
of CLRC, ARC, RITS, RDRC, CUL4-DDB 1 , HIRA, Clr6C, 
TRAMP, and SHREC, and in individual factors such as HP1 
homolog Swi6, NAD + histone deacetylase Sir2, CENP-B homolog 
Cbpl, and CHD1 remodeler Hrpl (Figure 1C). There were also a 
number of previously reported heterochromatin mutants that are 
listed in the library but were missed in our screen, such as agolA, 
rdplA, and dcrlA. We confirmed by PCR that the null mutation 
was missing in these strains, indicating that false negatives are 
more likely the result of incorrect deletions present in the library 
than due to methodological bias. 

Most interestingly, there were eight novel mutants identified in 
this screen, of which four are uncharacterized genes implicated in 
various steps of mRNA splicing: cwfl4 (SPBC24C6.il), dre4 
(SPAC13C5.02), cw/12 (SPBC32F12. 05c), and the human SRRM1 
homolog we named srml (SPCC825.05c). Each of these strains 
showed elevated levels of pericentric dh transcripts, a common 
phenotype of heterochromatin mutants, suggesting that these 
mutants indeed affected pericentric heterochromatin assembly 
(Figure ID). As cwfl4A showed the strongest phenotype among 
these four splicing mutants, we chose it as the focus of our 
subsequent experiments. 
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Figure 1. A genetic screen for nonessential genes required for pericentric heterochromatin silencing in fission yeast. (A) Schematic 
diagram of the otr::ade6 + -natMX6 reporter. (B) Workflow to introduce otr::ade6 + into the deletion library. A control screen was performed in parallel. 
The selection of only ade6-M21 0-hphMX6 progeny generates a uniform background for color development. White arrows denote a mutant that 
affects silencing, and yellows arrows denote a false positive that affects colony color independently of silencing. (C) List of mutants identified 
affecting pericentric silencing. Hit colonies were assigned color scores between 1 and 4, as indicated. (D) qRT-PCR analyses of transcripts derived from 
pericentric dh repeats, normalized to actl. Wild-type was set to 1. Error bars represent standard deviation of three experiments. 
doi:1 0.1 371 /journal.pgen.1 004334.g001 



Cwf14 is required for RNAi-mediated pericentric 
heterochromatin assembly 

We first confirmed by serial dilution analysis that cwfl4A cells 
containing atr::adeo formed white colonies similar to those of dcrlA 
(Figure 2 A). We also examined the effect of cwfl4A on silencing of 
another reporter inserted at the same location, otr::ura4 + [26]. The 



silencing of this reporter gene in wild-type cells allows them to 
grow on counterselective medium containing 5-fluoroorotic acid 
(FOA). Serial dilution analyses showed that cmfl4A has compa- 
rable silencing defects to dcrlA, as indicated by attenuated growth 
on FOA medium (Figure 2A). Moreover, chromatin immunopre- 
cipitation (ChIP) analyses showed increased enrichment of RNA 
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Figure 2. Cwf14 is required for pericentric heterochromatin assembly through the RNAi pathway. (A, E, F, and G) Serial dilution analyses 
to measure the expression of reporter genes and sensitivity to TBZ. (B and C) ChIP analyses of Pol II (Rpb1 ), H3K9me2 : and Swi6 levels at pericentric dh 
repeats and otr::ura4 + . Pol II ChIP was normalized to act1 gene, and H3K9me2 and Swi6 ChIP was normalized to acf? promoter. Error bars represent 
standard deviation of three experiments. (D) Northern blot analyses of siRNAs derived from pericentric dh repeats. 
doi:10.1371/journal.pgen.1004334.g002 



Polymerase II at pericentric regions in cwfl4A cells, indicating that 
the effect on pericentric transcript levels was at least in part due to 
increased transcription (Figure 2B). 

Further ChIP analyses showed strong reduction in levels of 
heterochromatin hallmarks such as H3K9me2 and Swi6 at 
otr::ura4*~ and to a lesser extent at the endogenous dh repeats in 
cwfl4A cells (Figure 2C). This pattern is similar to RNAi mutants 
such as dcrlA, but less pronounced than that of ck4A, suggesting 
that Cwfl4 might regulate RNAi-mediated heterochromatin 
assembly. Consistent with this idea, siRNAs derived from 
pericentric repeats were eliminated in cwfl4A cells, similar to 
dcrlA cells (Figure 2D). 



It was shown that a 1.6 kb fragment of pericentric dg 
repeats (termed L5) induces heterochromatin assembly and 
silencing of adjacent genes when inserted into an ectopic site in 
an RNAi-dependent manner [27,28]. Silencing is lost in 
cwfl4A cells as well, consistent with the idea that Gwfl4 is 
involved in RNAi-mediated heterochromatin assembly 
(Figure 2E). 

Since pericentric heterochromatin promotes proper loading of 
cohesin near centromeres to promote chromosome bi-orientation 
[29-31], most heterochromatin mutants show defects in chromo- 
some segregation [26] and sensitivity to the microtubule-destabi- 
lizing drug thiabendazole (TBZ) [32]. As expected, cw/14A was 
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also sensitive to TBZ (Figure 2A). Moreover, both the silencing 
and TBZ-sensitivity phenotypes were rescued by complementation 
of cwfl4A with a plasmid containing an intact copy of cwfl 4^ 
under the control of its endogenous regulatory elements 
(Figure 2F). Since cwfl4A cells had some growth defects, we also 
performed PCR-mediated random mutagenesis and isolated a 
mutant (cw/14-F26L) that resulted in silencing defects without 
significantly affecting growth (Figure 2G). 

To test whether cwfl 4 A affects H3K9 methyltransferase activity of 
CLRC, we used a system in which Clr4 is tethered via a fused Gal4 
DNA binding domain (GBD) to an ectopically integrated ade(f 
reporter adjacent to three copies of the Gal4 binding site (Sxgbs-adetf) 
[33] (Figure 3A). This tethering induces heterochromatin formation 
over a 6 kb locus, silencing transcription of the ade6 + reporter. RNAi 
factors are not required for diis silencing, consistent with the idea 
that RNAi is required for the targeting of CLRC to DNA repeats 
[33]. Dilution analysis showed that cwfl 4 A had no effect on Clr4- 
tethered silencing of 3xgbs-adeo ', and ChIP analyses showed that 
H3K9me2 was enriched at nearby loci 2 and 3 kilobases away in 
cwfl 4 A cells, similar to wild-type cells (Figure 3B). Collectively, these 
data suggest that Cwfl 4 is involved in heterochromatin formation at 
pericentric repeats through the RNAi pathway. 

Cwfl4 associates with the spliceosome 

Cwfl 4 is highly conserved across species. Its budding yeast 
homolog, Bud31p, co-purifies with the spliceosome [34], and 
BUD31A causes mis-splicing of the mRNAs of ARP2 and SRC1, 
two factors required for proper budding [35]. Cwfl 4 was initially 
identified in a purification of splicing factor Cdc5, together with a 
number of spliceosome components [36]. However, whether 
Cwfl 4 is a stable component of the spliceosome and involved in 
splicing has not been tested directly. In order to further specify the 
mechanism of Cwfl 4 action, we constructed C-terminalfy tagged 
versions of cwfl 4 at its endogenous locus. Cwfl4-GFP and Cwfl 4- 
myc were fully functional, as they did not show silencing defects of 
otr::adetf~ (Figure 4A). Imaging of Cwfl4-GFP showed that Cwfl4 
localizes predominantly in the nucleus (Figure 4B). Western blot 
analysis showed that Cwfl4-myc is a 40 kD protein. Moreover, 
the F26L mutation resulted in reduced Cwfl 4 protein levels, 
indicating that it is a partial loss-of-function allele (Figure 4C). We 
also performed immunoprecipitation of Cwfl4-myc followed by 
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MudPIT mass spectrometry to identify its interacting proteins. 
Most factors that co-immunoprecipitated with Cwfl 4, but not in a 
control purification, were components of the spliceosome, 
especially from subcomplexes NTC and U5 (Figure 4D and 
Table SI). This result corroborates data that Cwfl4 co-precipitates 
with Prpl7, Prpl9, and Cwf2, all members of U5-associated NTC 
[19,37], as well as a component of U5 snRNP, CwflO [38]. 

Cwf14 is involved in splicing of specific mRNAs 

Previous work suggests that the spliceosome is involved in 
heterochromatin assembly through tethering RDRC to pericentric 
transcripts [15]. This is because RDRC component Cidl2 
associates with the spliceosome [15,20], and no splicing defects 
were observed in the well-characterized splicing substrate tbpl in 
splicing factor mutants when silencing defects were apparent 
[15,17]. However, purifications of spliceosome components, 
including our purification of Cwfl 4, have not identified any other 
heterochromatin components [36,37,39] (Table SI). These results 
indicate either that the physical connection between the spliceo- 
some and RDRC is very weak, or that Cidl2 is present in two 
separate complexes, RDRC and spliceosome. It remains a 
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Figure 4. Cwfl 4 associates with the spliceosome. (A) Tagged 
versions of Cwfl 4 are functional as indicated by serial dilution analysis. 
(B) Imaging of cells expressing Cwf14-GFP, with DAPI to stain nucleus. 
Scale bar, 10 ^m. (C) Western blot analysis of Cwf14-myc expression. (D) 
Mass spectrometry analyses of purified Cwf14-myc complex. Factors 
were binned into annotated subcomplexes. 
doi:10.1371/journal.pgen.1004334.g004 
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possibility that splicing factors regulate the correct processing of 
mRNAs involved in RNAi-mediated heterochromatin assembly. 

In order to test whether cwfl 4 A has a general splicing defect, we 
performed RNA-seq analyses of total cellular RNAs from wild- 
type, dcrlA, and cwfl4A cells. The gene expression profiles of 
cwfl4A and dcrlA showed strong overlap of significantly up- 
regulated genes (Figure 5A and Tables S2 and S3), consistent with 
the idea that Cwfl 4 functions in the RNAi pathway. We found 
that cwfl 4 A indeed resulted in intron retention of a portion of 
genes (Figure 5B and Table S4), consistent with the finding that it 
is associated with the spliceosome. The majority of introns were 
properly processed, which indicates that Cwfl 4 only moderately 
affected the activity of the spliceosome. Interestingly, many RNAi 
factors contain introns. Although unbiased ranking of exon-exon 
junction ratios between wild-type and cwfl 4 A RNAs did not show 
preferential enrichment of RNAi or heterochromatin assembly 



factors (Table S4), significant intron retention was observed within 
mRNAs from agol, arb2, ersl, and dshl in cwfl4A cells as compared to 
those in wild-type cells, despite similar sequencing depths genome- 
wide and similar levels of each gene transcript in both samples 
(Figure 5C). RT-PCR analyses confirmed that these mRNAs were 
indeed spliced inefficiently in cwfl 4 A as well as in cwfl 0-1 and prplO-1 
cells (Figure 5D). Moreover, Western blot analysis showed strong 
reduction of protein levels of Flag-Agol, the RNAi factor most 
severely affected by cwfl 4 A, in both cwfl 4 A and cwfl 0-1 strains 
(Figure SI), indicating that the mis-splicing of agol mRNA, and 
possibly other RNAi factors as well, resulted in altered protein levels. 

Previously, it was shown that introducing cDNA versions of 
agol or hrrl + failed to rescue silencing defects oiprplO-1 [15]. We 
also generated a cDNA version of agol at its endogenous 
chromosomal location and under the control of its endogenous 
regulatory elements {agol ::cDNA). This cDNA construct showed 
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no defects in silencing of otr::uraf, indicating that this replacement 
created a functional Agol protein (Figure S2). However, 
agol + ::cDNA was unable to rescue otr.mraV silencing defects and 
TBZ sensitivity of cwfl4A (Figure 6A), even though it restored 
Flag- Agol protein levels (Figure SI). We reasoned that the 
inability of agol + ::cDNA to rescue avfl4A defects is because other 
RNAi factors are also mis-spliced. We thus generated cDNA 
versions of arb2* and ersl + at their endogenous chromosomal loci, 
which were both functional (Figure S2). Neither arb2 + :xDJVA nor 
ersl + ::cDNA alone showed any effect on silencing of otr::ura-f~ in 
avfl4A cells (Figure 6A). However, when combinations of two 
cDNAs were introduced together into cw/14A cells, there was a 
detectable rescue of silencing defects and TBZ sensitivity, and the 
effect was stronger when all three cDNAs were introduced 
(Figure 6A). Further ChIP analysis showed that both H3K9me 
and Swi6 protein levels were significantly increased at both 
otr::ura4 + and pericentric dh repeats in cwfl4A cells supplemented 
with agot, arb2 + , and erst cDNAs (3cDNAs) (Figure 6B). 

We also found that introducing agol + ::cDNA was sufficient to 
rescue the silencing defects of otrr.adeo associated with cwflO-1 
(Figure S3A). Moreover, there is a significant increase in both 



H3K9me and Swi6 levels at both otry.adef? and dh repeats in cwflO- 
1 agol ::cDNA cells (Figure S3B). Altogether, these results suggest 
that mutations in different splicing factors affect the splicing of 
diverse RNAi factors to regulate heterochromatin assembly at 
pericentric regions. 

Thus our results clearly demonstrated that splicing factors 
mainly exert their effects on pericentric heterochromatin 
assembly by promoting the proper splicing of RNAi factors. 
However, we caution that the rescue of heterochromatin 
silencing defects of cwfl4A cells was incomplete even with 
agol , arb2 + , and ersl + cDNAs. This probably reflects the 
requirement of Cwfl4 for the proper splicing of other RNAi 
factors such as dshl + , rdpl + , arbl + , hrrl + , or some unidentified 
factors involved in heterochromatin assembly. Such an idea is 
supported by the fact that the rescue of cwfl 4 A silencing defects 
correlated with the number of cDNA constructs that were 
introduced. It is also possible that splicing factors have a direct 
contribution in pericentric heterochromatin assembly. However, 
the strong rescue of pericentric silencing in cw/14A cells with cDNA 
constructs suggests that the regulation of RNAi factor splicing is a 
major role of splicing factors in this process. 
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Cwfl4 regulates the proper splicing of tpzl to control 
telomere silencing 

We also found that telomere shelterin component tpzl + showed 
a strong reduction in exon-exon junction sequencing reads in 
cwfl4A cells relative to wild-type cells, and this phenotype was 
confirmed by RT-PCR analyses in cwflO-1 and prpl()-l as well 
(Figure 7A). A C-terminally Flag-tagged version of Tpzl [40] 
affected silencing of a reporter gene inserted near telomere repeats 
(Figure S4), suggesting that Tpzl is required for telomere silencing. 
We then analyzed the effect of splicing factor mutations on 
silencing of a reporter inserted near telomere repeats of 
chromosome two (Tel2::urat) [41]. Indeed, cwfl4A resulted in 
silencing defects at this reporter gene (Figure 7B), accompanied by 
a reduction of H3K9me levels at tlhl, which is embedded at 
telomeric heterochromatin, as well as the accumulation of tlhl 
transcripts (Figure 7C and 7D). Both cwflO-1 and prplO-1 cells 
showed increased tlhl transcript levels, indicating that loss of 
telomeric silencing is a general feature of splicing mutants 
(Figure 7C and 7D). Because multiple DNA sequences contribute 
to heterochromatin assembly at telomeres, including tlht, 
telomere associated sequences (TAS), and terminal TEL repeats 
[8], we further tested silencing at TEL::ade6* ', which is inserted on 
the mini-chromosome Chi 6 adjacent to telomere repeats [42]. We 
found that cwfl4A, cwJ14-F26L, and cwflO-1 resulted in loss of 
silencing of this reporter (Figure 7E). Interestingly, replacement of 
tpzt with its cDNA partially rescued telomeric silencing phenotypes 
of the cwfl 4-F26L and cwfl 0-1 mutants (Figure 7E and 7F). We 
could only marginally rescue the telomere silencing defects of 
cwfl 4A cells with tpzt '-'-cDNA (Figure 7F and data not shown). Thus 
it is possible that additional factors involved in telomere silencing 
might also contain introns and depend on the spliceosome to 
properly regulate their splicing. Alternatively, splicing factors might 
affect telomere silencing through additional mechanisms. Nonethe- 
less, these results demonstrate that inefficient splicing of tpzl 
contributes to telomeric silencing defects in splicing mutants. 

Discussion 

The formation of heterochromatin requires RNAi-mediated 
processing of repeat-derived transcripts and the targeting of 
histone modifying activities to repeat regions, leading to H3K9me 
and the recruitment of HP1 proteins. It has been shown that 
splicing factors are required for RNAi-mediated heterochromatin 
assembly in fission yeast, although the mechanism by which they 
are involved is not well characterized. The previously prevailing 
model was that the spliceosome physically associates with RNAi 
factors to regulate heterochromatin assembly, rather than acting 
through its splicing activity [15,16]. One of the main evidences for 
this idea is that splicing mutants show severe silencing defects even 
though the splicing of thpl mRNA was not affected [15,17]. 
However, whether these splicing factor mutants selectively affect 
the splicing of RNAi factors has not been rigorously tested. 

Our RNA-seq analyses showed prominent intron retention of a 
subgroup of mRNAs in cwfl 4 A cells, even though the majority of 
introns (including those of tbpl) are still properly processed 
(Figure 5B and Table S4). Interestingly, a number of key RNAi 
factors were among the list of strongly mis-spliced introns, a result 
that is further corroborated by RT-PCR analyses of a selective set 
of RNAi factor mRNAs in other spliceosome mutants such as 
cwflO-1 and prplO-1 (Figure 5D). Most importandy, we found that 
introducing a combination of cDNAs of RNAi factors significandy 
alleviated pericentric heterochromatin defects associated with 
cwfl 4 A and cwfl 0-1 (Figure 6 and S3). Thus splicing factors 
regulate the proper splicing of RNAi factors, which is a major, if 



not sole, contributor to heterochromatin assembly defects in 
splicing mutants. That introducing tpzl* "r.cDNA was able to 
partially rescue telomere silencing defects associated with splicing 
factors further supports the idea that mis-splicing of heterochro- 
matin factors is the reason splicing factor mutants show 
heterochromatin assembly defects (Figure 7E and 7F). 

It is noteworthy that cwfl 4 A has only moderate splicing defects, 
with some introns show very strong sensitivity, whereas most 
others show littie to no defects (Figure 5B and Table S4). This 
raises the question of whether specific intron sensitivity is a result 
of introns that are inherently difficult to splice. Consistent with this 
idea, our RT-PCR analyses showed prominent unspliced precur- 
sor mRNAs of RNAi factors even in wild-type cells (Figure 5D). It 
is also a striking pattern that heterochromatin factors in the same 
complexes tend to either have or not have introns. For example, 
many members of RNAi, such as agol + , arbt, arb2f, ersl + , dshl + , 
rdpl + , and krrt ', have introns, but none of CLRC [elm, rikl + , 
raft, raft, cult, stct), SHREC (clrt, clr2 + , clr3 + , mitt), or swift 
have any introns, raising the possibility of selective regulation of 
the RNAi pathway through general or specific changes in splicing 
efficiency. In fact, an analysis of splicing changes during meiosis 
showed that one intron of arbt is induced to be spliced, while a 
different intron exhibits splicing repression [43]. Since splicing 
factors have been identified in screens that affect RNAi-based 
processes in worms, flies, and plants [44-48] , it seems a conserved 
mechanism that proper splicing of the mRNAs of RNAi factors 
regulates the efficiency of RNAi inside the cell. 

Materials and Methods 

Genetic screens of the fission yeast deletion library 

The Bioneer library is composed of strains of mixed endogenous 
ade6 alleles: ade6-M216, which forms pink colonies on low adenine 
medium, and ade6-M210, which give rise to red colonies, similar to 
adedA. To avoid the complication of ade6-M216 alleles in our screen, 
we included in the query strain an ade6-M2 1 0-mCherry allele linked to a 
hphMX6 cassette, which confers resistance to the antibiotic hygro- 
mycin B, allowing us to generate a uniform ade6-M210 background. 
Screens were carried out according to previous protocols [49], with 
slight modifications. Yeast strains arrayed in 384 strains/ plate format 
were pinned on YES agar medium containing 100 ug/ml G418, and 
otr::ade(f-natMX6 ade6-210-hphMX6 and ade6-210-hphMX6 strains 
were pinned on YES+100 ug/ml hygromycin B. After two days 
growth, strains were mated on SPA agar medium. Plates were 
incubated at 25°C for 3 days, then 42°C for 3 more days to kill 
vegetative cells. Strains were then germinated and correct genotype 
selected by pinning to YES+GNH (G418, Nat, and Hygromycin) or 
YES+GH (G4 1 8 and Hygromycin) medium, and then pinned to YE 
(low ade) medium for color readout. 

Fission yeast strains and genetic analyses 

Cwfl4-GFP and Cwfl4-13myc strains were constructed by a 
PCR-based module method, cwfl 4 A, cwfl 2 A, dre4A, and srmlA were 
derived from the Bioneer fission yeast deletion library, verified via 
PCR, and backcrossed. The p-ewflt plasmid was constructed by 
insertion of a PCR product containing cwfl t promoter and coding 
region into SphI and Xmal sites of pREP41. The cwfl 4-F26L-myc 
mutant strain was constructed by integrating a cwfl f ' -myc-kanMX6 
cassette amplified by error-prone PCR (high dNTP concentration) 
into the endogenous cwflt locus. agotr.cDMA, arb2f::cDNA, and 
erstr.cDMA strains were constructed by replacing the endogenous 
genes with intronless cDNA versions. All were sequenced to confirm 
full replacement and lack of mutation, cwfl 0-1 and prplO-1 strains 
were a kind gift from Robin Allshire. Genetic crosses were used to 
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Figure 7. Splicing factors regulate the proper splicing of shelterin component Tpzl. (A) RT-PCR analyses of RNA with primers flanking 
introns. (B, E) Serial dilution analyses to measure reporter gene expression. (C) ChIP analyses of H3K9me2 levels, normalized to act1 promoter. Error 
bars represent standard deviation of three experiments. (D, F) qRT-PCR analyses of transcripts derived from tlhl, normalized to act1 gene. Wild type 
was set to 1. Error bars represent standard deviation of three experiments. 
doi:1 0.1 371 /journal.pgen.1 004334.g007 



construct all other strains. The genotype of strains is provided in 
Table S5. For serial dilution plating assays, ten-fold dilutions of 
mid-log-phase culture were plated on the indicated medium and 
grown for 3 days at 30°C unless otherwise indicated. 

RNA analyses 

Total cellular RNA was isolated from log-phase cells using 
MasterPure yeast RNA purification kit (Epicentre) according to 
the manufacturer's protocol. Quantification with qRT-PCR was 



performed with Power SYBR Green RNA-to-CT one-step Kit 
(Applied Biosystems). RNA serial dilutions were used as templates 
to generate the standard curve of amplification for each pair of 
primers, and the relative concentration of target sequence was 
calculated accordingly. An actl fragment served as reference to 
normalize the concentration of samples. Sequence of DNA oligos 
is provided in Table S6. For RNA-seq, purified RNA was 
prepared by TruSeq Stranded Total RNA Kit (Illumina), which 
includes rRNA depletion and chemical fragmentation. Index 
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adapters were added to allow for multiplexing. Paired-end 
sequencing with 100 bp read lengths was performed on Illumina 
HiSeq. Mapping was performed with the Tuxedo Suite consisting 
of Bowtie, TopHat, and Cufflinks. For cwfl4A, 57,521,513 reads 
were obtained, and 83% mapped to the genome via Bowtie. For 
WT, 47,636,518 reads were obtained, and 82% mapped. For 
dcrlA, 64,972,007 reads were obtained, and 83% mapped. RNA- 
seq data have been deposited to the Sequence Read Archive 
(http://www.ncbi.nlm.nih.gov/sra/) with accession number 
SRP040479. For the dot plot, exon-exon junction ratios were 
filtered to remove several types: 1) junctions which mapped zero 
times in any sample (possible mapping noise), 2) junctions whose sum 
in WT, cwfl 4 A, and dcrlA was less than 30 (to avoid randomness due 
to small sample sizes), and 3) ratios whose values were greater than 
150 (to focus the diagram on splicing reduction). Northern blot of 
siRNAs was performed as described previously [50]. 

Chromatin immunoprecipitation (ChIP) analysis 

ChIP analyses were performed as described previously [51]. 
Antibodies used were H3K9me2 (Abeam 1220), Swi6 [52], and 
RNA Pol II (Covance 8WG16). Quantification with qPCR was 
performed with Maxima SYBR Green/ROX qPCR Master Mix 
(Thermo). Enrichment was calculated as: relative levels in ChIP/ 
relative levels in total DNA. An actl promoter fragment was used 
as a control for normalization unless otherwise indicated. 
Sequence of DNA oligos was provided in Table S6. 

Protein purification and mass spectrometry analysis 

Exponentially growing yeast cells were harvested, washed with 
2xHC buffer (300 mM HEPES-KOH pH 7.6, 2 mil EDTA, 
100 mM KC1, 20% glycerol, 0.1% NP-40, 2 mM DTT, and 
protease inhibitor cocktail (Roche)) and frozen in liquid nitrogen. 
Crude cell extracts were prepared by vigorously blending frozen 
yeast cells with dry ice using a household blender, followed by 
sonication and incubation with 30 ml lxHC buffer containing 
250 mM KC1 for 30 min. The lysate was cleared by centrifugation 
at 82,700 xg for 3 hours. The supernatants were incubated with 
50 JU.1 of C-myc antibody (Sigma C3956) overnight, for three hours 
the next day with 50 |il protein G agarose beads, washed eight 
times with lxHC containing 250 mM KC1, then two times with 
the same buffer without NP-40. For mass spectrometry analysis, 
bound proteins were eluted with 2 x 100 (U.1 of 50 mM Tris pH 7.5, 
5% SDS, 5% glycerol, 50 mM DTT. MudPIT mass spectrometry 
analysis was performed as described previously [53]. 

Supporting Information 

Figure SI Western blot analysis of Flag-Ago 1 protein levels. Cell 
lysates were first immunoprecipitated with Flag-agarose beads and 
Western blot analyses were performed with a Flag antibody. 
(PDF) 
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